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Abstract. We study systems with two symmetric absorbing states, such as the 
voter model and variations of it, which have been broadly used as minimal neutral 
models in genetics, population ecology, sociology, etc. We analyze the effects of a key 
ingredient ineluctably present in most real applications: random-field-likc quenched 
disorder. In accord with simulations and previous findings, coexistence between the two 
Q competing states/opinions turns out to be strongly favored by disorder in the standard 

voter model; actually, a disorder-induced phase transition is generated for any finite 
-H system size in the presence of an arbitrary small spontaneous-inversion rate (preventing 

absorbing states from being stable). For non-linear versions of the voter model a 
general theory (by AlHammal et al.) explains that the spontaneous breaking of the 
^^ up/down symmetry and an absorbing state phase transition can occur either together 

^^ or separately, giving raise to two different scenarios. Here, we show that he presence 

■^- of quenched disorder in non-linear voter models does not allow the separation of the 

up-down (Ising-like) symmetry breaking from the active-to-absorbing phase transition 

in low-dimensional systems: both phenomena can occur only simultaneously, as a 

• • consequence of the well-known Imry-Ma argument generalized to these non-equilibrium 

. J_< problems. When the two phenomena occur at unison, resulting into a genuinely non- 

j^ equilibrium ("Generalized Voter") transition, the Imry-Ma argument is violated and 

3 the symmetry can be spontaneously broken even in low dimensions. 
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1. Introduction 

In Nature, thermodynamic equilibrium is the exception rather than the rule. Almost 
the entire world of observable collective phenomena -from the formation of galaxies and 
stars to the self organization of communities of insects or neurons in the brain- can be 
best adressed within an out-of-equilibrium framework [JJ. Owing to this, huge efforts 
have been devoted to the study of the statistical mechanics of nonequilibrium systems 
and their applications in physics, biology, ecology etc. Understanding the fundamental 
laws governing the statistics and the dynamics of these systems is a formidably complex 
task; still, a great variety of them -even if a priori very different in their nature- show 
similar underlying statistical features that can be captured by rather simple probabilistic 
toy models, a few of which have become paradigmatic [21 H3 HI E] • 

The Voter Model (VM) [2] occupies a preponderant position in this context. It was 
independently proposed in various research fields to study neutral genetic drift in an ideal 
population [51 E], the competition for territory between two countries [8], spreading of 
infectious diseases [9] or language competition. It is one of the few known interacting- 
particle models which is exactly solvable in any spatial dimension [21 [JU] and is now 
extensively employed and studied -even if with different names- in population genetics 
[IT] , ecology [121 [13], sociology [H], and linguistics [HI [16] among other disciplines. The 
VM is simply defined: In a community of individuals (voters) arranged on the vertices 
of a regular lattice (or, more generally, upon an arbitrary network [13]) each voter 
has one out of two equivalent states (opinions) -generalizations to more states being 
straightforward- and at every discrete time step the dynamics is defined as follows: 

(i) a voter is randomly selected, 

(ii) its opinion is substituted by the one of a randomly chosen nearest neighbor, 

(iii) the process is iterated ad infinitum or until consensus (all voters sharing the same 
opinion) is reached. 

Some relevant features which determine the global behavior of the model and its 
observables are [21 CH El H] : 

• the dynamics is completely neutral, that is, the opinions' labels can be switched 
leaving the dynamics unaltered, this makes the VM the most basic model to analyze 
neutral theories in population genetics (where the states correspond to "alleles") 
and in ecology (in which the states are different "species"); 

• it is a linear model in the sense that the probability for a single voter to change 
its opinion, that is, the flipping probability f, increases linearly with the number 
of discordant neighbor voters. This allows for an exact moment-closure of the 
dynamical equations for the n— point correlators; 

• it has no free parameters and since it lacks of any type of characteristic (length or 
time) scale it seats, by definition, at a critical point; 

• it is characterized by a purely noise-driven diffusive dynamics without surface 
tension at the boundaries between different domains of opinion; 
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• it exhibits two Z2-symmetric absorbing states, from which the system cannot 
escape. 

Some key properties can be summarized as follows [TH]: For lattices of dimension 
D < 2 the system reaches one of its two absorbing states with probability one in the 
thermodynamic limit; instead, for D > 2 the active state -with two coexisting opinions 
- lasts indefinitely. For finite-size systems consensus is eventually reached almost surely, 
i.e. with probability one, since it is the only state in which the dynamics ceases. If N 
is the total number of voters, the mean time to reach consensus, T/v, scales as N 2 in 
D = 1, iVlogiV in D = 2, and as iV for D > 3. Similarly, in the infinite-size limit, the 
2-point correlation function G(r, t), measuring the probability that two voter at distance 
r at time t have the same opinion, approaches 1 at arbitrary fixed distance r at large 
times for D < 2, while it goes G(r, t) ~ r 2 ~ D for D > 2, consistent with the lack of 
consensus. 

Nonlinear generalizations of the VM have been considered, for example, in neutral 
ecology, where the VM is the simplest model mimicking stochastic species competition. 
Nonlinearity in this case stems from a dependence of the dynamical rules on the density 
of individuals of each specific species in a local neighborhood [201 EB [22J [23] . Actually, it 
is well understood that negative density-dependence (implying that a locally infrequent 
species has a competitive advantage) significantly favors species coexistence. Similarly, 
a set of models labeled "with heterozigosity selection" where introduced in population 
genetics (having locally different coexisting "alleles" is favored or, in other words, 
minorities tend to be preserved) [21]. Different non-linear variations of the VM have been 
reported to exhibit a novel type of phase transition, lying in the so-called Generalized 
Voter (GV) universality class (see [171 [18], [251 I2H1 [27] and refs. therein). This class is 
characterized by a critical point separating and active from an absorbing region, whose 
scaling right at criticality coincides with that of the original VM. Moreover, its dynamics 
is dual to a family of models of branching and annihilating random walkers with parity 
number conservation [TOl [28] . 

A key ingredient likely to strongly influence the dynamics of real systems is quenched 
disorder. Disorder -which is unavoidable in Nature- can be defined naively as an 
intrinsic component of randomness in the interaction among the "micro-constituents" 
or in the topological structure over which the dynamics takes place. It is well known 
from statistical mechanics that quenched disorder can have a dramatic impact in 
equilibrium systems [291 HSU EH [32]. Recently, there has been a growing interest in 
studying the effects of disorder in genuinely nonequilibrium models and in particular 
in models with absorbing states [331 EU ESI [36]. Frachebour, Krapivsky, and Redner 
[37] studied the influence of quenched disorder in the form of impurities for a model 
of catalysis with two symmetric absorbing states, showing that a non-trivial steady 
state emerge. More recently, Masuda et al. [3B1 EH] showed that quenched (random- 
field like) disorder -creating an intrinsic preference of each individual for a particular 
state/opinion- hinders the formation of consensus, hence favoring coexistence. Actually, 
the presence of just a few different "zealots" -not allowed to change their intrinsic 
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state- suffices to prevent consensus [40]. Along similar lines, Pigolotti and Cencini 
[41] analyzed in the context of neutral ecology a version of the VM in which at each 
location there is an intrinsic preference for one particular species, leading to mixed states 
(no consensus/monodominance) lasting for times that grow exponentially with system 
size. By studying a similar model, Barghathi and Vojta [42] have very recently stressed 
that -contrarily to what happens in equilibrium systems, where a well known (Imry- 
Ma) argument precludes symmetries to be spontaneously broken in low-dimensional 
systems in the presence of quenched random fields [43J, intrinsically non-equilibrium 
phase transitions, such as those in the GV class, do persist in low-dimensional systems 
{D = 1) in the presence of random fields, even if with a different type of critical behavior 
[42] . Despite of all these results, a complete and coherent theoretical framework to 
understand the effects of disorder in VM-like systems is still missing. 

In this paper, aimed at shedding further light on this problem, we consider a 
VM where each voter experiences an intrinsic tendency to align with a particular 
opinion. We propose a systematic study of this type of models with an approximate 
analytical description of the underlying stochastic dynamics complemented with 
computer simulations. 

The paper is organized as follows: In the next section we introduce a model with two 
symmetric absorbing states and quenched random fields. In section 3 we first discuss the 
mean-field (deterministic) limit of the model and obtain an approximate solution taking 
into account finite-size (stochastic) effects. We find that disorder keeps the system 
away from the absorbing boundaries and -in agreement with previous findings- the 
typical time to fixation is exponential in the system size instead of linear as in the VM. 
Section 4 discusses the connections with the GV universality class. Section 5 extends 
the analysis to a class of nonlinear voter model, while in Section 6 concluding remarks 
and considerations are presented. 

2. The model 

2.1. Definition of the model and notations 

We consider a Voter Model defined on a D-dimensional lattice (A C Z D ) with N = L D 
sites, denoted by i,j,... [21 OH]- At each site i resides a binary or spin variable 
o~i G {+1,-1} and a random binary field Tj to which Oi is locally coupled, favoring its 
alignment with the field. The values of r are quenched, that is, a particular realization 
of the disorder is extracted and does not change during the dynamics. To ensure the 
global up-down (plus/minus) symmetry we partition the lattice into two disjoint sets of 
the same size, namely A = A + U A~, with A 1 * 1 = {i G A : r» = ±1} and |A + | = |A~|. 
In finite spatial dimensions, for large enough systems, this constraint can be relaxed 
to be satisfied just on average over the lattice: we take tj as i.i.d. random variables 
taking values in {+1,-1} with uniform probability. At each site the coupling-strength 
between the spin and the random field is controlled by the free parameter e G [0, 1] , 
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Figure 1. Cartoon of the microscopic dynamical rules of the model. (Left) The 
underlying lattice is represented by black/white nodes representing the quenched 
random field pointing upward/downward. A central "spin", i, with associated random 
field pointing upward is removed and replaced by (Right) a nearest-neighboor state; 
the relative transition rates to chose any (only two of them are depicted with red 
dashed arrows) depends on whether it is aligned or not with tj: 1 + e if it becomes 
aligned, 1 — e if it does not. 



where e = stands for the uncoupled (pure VM) case and e = 1 implies that each spin 
remains frozen in the direction of its random field. The model is completely defined -in 
its continuous-time version- by specifying the transition rates W for a generic spin i, 
see Figure [TJ namely: 

W{a % -► -a,) = ^p £(1 - aw), (1) 

jedi 

where di is the set of nearest neighbors of % and z is the lattice coordination number 
(z = 2D, for a regular square lattice). The term independent of e describes the standard 
VM dynamics [TD], while the second term is proportional to e and the flipping probability 
is enhanced or reduced depending on whether the spin becomes aligned or not with its 
random field. 



2.2. Mapping onto a birth-death Fokker-Planck equation 

As a first step to construct a mean-field solution, let us consider the dynamics on 
a complete graph, that is, each spin connected to every other spin (in eq. ([IJ this 
corresponds to di = {j ^ i}). The macroscopic state of the system is univocally 
determined by the value of two variables, x and y, which represent the fraction of up 
and down spins aligned with their corresponding random fields, respectively: 



N 



i=l 

N 



T; 



T; 



+ 1) 



(2) 



d t P(x,y,t) = d x 



-A X P + ^d x (B x P) 



+9y 



~AP + 7^-MByP) 



(4) 
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These two variables are defined in the interval [0, 1/2] and, since the total number N of 
spins is constant, the total fraction of up and down spins (X and Y respectively) are 
readily obtained as X = 1 — Y = 1/2 + 1-1/6 [0,1]. The global magnetization is given 
by (ft = 2(x — y). We can map this spin model onto a birth-death process considering the 
Master Equation (ME) for P(x, y, t), the joint probability of having at time t a fraction 
x and y of up and down spins aligned with their local field, respectively. In this (fully 
connected or mean-field) version of the model, P(x,y,t) evolves through discrete steps 
x — > x' = x±l/N and y —y y' — y±l/N, with transition rates given by (see [38 | 1391 I^T] ) 

W b x = W{x -> x + l/N) = {l + e){\-x){\ + x-y) 
W* = W(x ->x- l/N) = (1 - e)x{\ -x + y) 
W b y = W{y-±y+ l/N) = (1 + e){\ - y){\ + y - x) 
W d y =W{y^y- l/N) = (1 - e)y{\ -y + x), 

The standard Kramers-Moyal expansion [H] leads to the Fokker-Planck approximation 
of the original ME for the evolution in time of P(x, y, t), that we write as 

1 

2N~ 

where Ai = Wf — W\ represent the drift terms and B>i = Wf + W\ the diffusion terms, 
% — x, y, and time has been rescaled in units of l/N. 

3. Steady state Analysis 

3.1. Deterministic limit 

For the time being we focus only on the limit iV — y 00, when the diffusion terms can 
be safely set to zero and the dynamics becomes deterministic, namely x = A x and 
y = A y . We perform a change of variables [38l [39] that will be useful for the later 
analysis by defining E = x + y G [0, 1] and A = x — y G [—1/2, 1/2], from which the 
global magnetization can be written as <fi = 2 A. In this notation, the deterministic 
equations become 

A = eA(l - 2E) 

(5) 
E= i(l + e)-S-2eA 2 . 

The analysis of the dynamical system described in Eq. (fsj) gives already some interesting 
results (as already outlined in [3S1 ESI SI]). For e = one obtains a line of stable fixed 
points at S = 1/2 (and arbitrary A), recovering the VM results and hence the system 
always reaches an absorbing state when fluctuations are considered, also in the infinite- 
size system. Instead, for e > 0, the phase portrait changes dramatically and the line 
of fixed points breaks into three fixed points: two of them are unstable corresponding 
to the absorbing states of the VM dynamics, at E = 1/2 and A = ±1/2, while the 
third, at S = 1/2(1 + e) and A = 0, is stable and corresponds to an active state with 
zero magnetization, i.e. a phase of symmetric coexistence of the two opinions. Thus 
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the infinite-size limit for this latter case is always in the active phase, but what are the 
effects of fluctuations when the size is finite? We will answer this question in the next 
sections. 

3.2. Finite-N approximate solution and role of the stochastic noise 

For finite-size systems fluctuations cannot be neglected and the system is expected to 
fluctuate around the deterministic stable fixed point. We expect a priori that only a 
large collective deviation can bring the system to one of the absorbing states which 
is, nevertheless, ineluctably reached. To study this we consider the Fokker-Planck eq. 
Q. It is easy to verify -by computing cross-derivatives- that this equation does not 
admit a potential solution for the stationary probability distribution jH]. The lack of 
a stationary potential reflects the intrinsically non-equilibrium nature of the problem. 
A possible strategy would be to construct non-differentiable non-equilibrium potentials 
following the strategy in jl5l H6] . Instead, here, we follow a simpler solution by seeking 
for a suitable (adiabatic) approximation allowing us to reduce the problem to a one- 
variable one |%7] . 

Let us consider e C 1, then Eq. (p)]) have two different characteristic relaxation 
times: £ relaxes in a time 0(1) whereas A in a much longer time-scale (9(e _1 ). Thus, 
one can assume that the system first relaxes to the nullcline orbit £ = and then the 
dynamics is constrained to take place uniquely on such a one-dimensional manifold. 
Consequently, we adopt a quasi-steady-state scheme where the variable E is substituted 
by its value in the nullcline orbit, namely 

£^£ = ^(l + e)-2eA 2 (6) 

and it is treated like a deterministic quantity, that is, fluctuations in its direction are 
discarded. 

Within this approximation the Fokker-Planck equation for the probability 
distribution of A, or equivalently <ft, V(<f),t), is obtained with a change of variables in 
Eq.Q from (x,y) to (A, E) and substituting the variable £ in the remaining equation 
with S of Eq. (pL The diffusion term for the variable S is neglected and we are left 
with the following 1-dimensional FP equation 

V(<f>, t) = -fy [A(<j>)V(<j>, t)} + \d\ [B{t/>)V(<f>, t)] , (7) 



2 



with 



A{4>) = -%0(l - <P 2 ) 

(8) 
S(0) = i(l-e 2 )(l-0 2 ). 

Eq. ([7]) is equivalent to the following Langevin equation in the Ito prescription [H] 

= A(<f>) + yfjmv{t) 

, (9) 

= -40(l-0 2 ) + ,A(l- e 2)(l-02 )r?(t)) 
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where rj is a (^-correlated in time gaussian white noise with zero mean. Let us emphasize 
that the main effect of the quenched disorder is to generate a deterministic force which 
stabilizes the opinion-coexistence state, = 0. In the limit e — > we recover the 
purely noise-driven VM dynamics [17J, while in the opposite limit e — > 1, the dynamics 
is purely deterministic and the spins align with their corresponding random fields (i.e. 
<p = on average). For < e < 1 the stationary solution V s (4>] e) is formally given by 
the zero-current condition 



from which 



J 



Vs 



-A(<f>)P(M + -d 4 ,[B(<l>)V(<i>,t)] 







1 



o ft> d r A{x) 
_ e z J ax B{x) 



(10) 



;n) 



' ZB(<j>) 

Z is supposed to be the normalization constant of V s {4>\ e), but since the diffusion term 
B(<p) — >■ when <p — > ±1 and the exponential stays finite, the probability distribution is 
not normalizable (as corresponds to the probability distribution collapsing to one of the 
absorbing states [18J). Therefore, as expected, for any finite value of N the only steady 
state is an absorbing/consensus one: coexistence is always killed on the large time limit. 



V(6») 
2.0 




Figure 2. Potential V(9) of Eq. |l4| ) corresponding to three values of e at fixed 
N = 200: e = 0.05 (lower blue curve) corresponds to the absorbing phase, e = 0.15 
(red central curve) and e = 0.2 (upper green curve) correspond to two examples of 
intermediate phase. The black curve corresponds to the parameters e = 0.2, iV = 250. 
Note how the potential well becomes deeper as e or N increase, diverging in the 
e > 0, N -> oo limit. 



The understand this problem even more transparently, we perform a change of 
variables on the Fokker-Planck equation ([7]) such that its corresponding Langevin 
equation ([9]), characterized by a state-dependent (multiplicative) noise, becomes a new 
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state-independent equation, i.e. the noise becomes additive rather than multiplicative. 
A suitable transformation is [49] 150] 

9 = — arcsin0 (12) 

a 



with a = y^f (1 — e 2 ), which leads to the following Langevin equation 
a dV(9) , t 



where 



d9 s ' 

1 e 2 



(13) 



V{9) = ^ log(cos(a#)) + — sm 2 {a6). (14) 

2 4cr 

is a potential and £ is a standard 5-correlated Gaussian white noise. Observe that 

the first term destabilizes the = (coexistence) solution, while the second, disorder- 



induced, one stabilizes it (indeed, the second term in eq. ( 14 ) can be expanded around 
9 = 0, leading to a parabolic potential around the origin; i.e. disorder creates an effective 
potential whose minimum corresponds to the opinion-coexistence state). V(9) is shown 
in Figure [2] for some particular values of the parameters e and N. We can see that there 
exists a critical value of e, e c (to be computed later) below which the potential effectively 
pushes the system toward the absorbing boundaries. Going above this critical value a 
local minimum appears in the configuration of zero magnetization = 0. We expect 
then that, for e > e c the time needed to reach the absorbing state will be exponential 
in the height of the potential barrier, due to the Arrhenius law [44]. 

As e and/or N increase, the basin of attraction of the minimum at <fi = 
becomes larger and deeper, that is, the two symmetric maxima of the potential become 
closer to the corresponding absorbing states and the barriers become higher, in such 
a way that the time needed to escape the barrier becomes much longer (actually, long 
enough to make it unaccessible to computer simulations). We can then identify three 
different regimes that we call the absorbing, intermediate (quasi-active) and active phase, 
respectively. In the absorbing phase symmetry is broken and one of the two possible 
states of consensus is reached with certainty; the active phase is characterized by a 
coexistence of both states (it survives to fluctuations only in the infinite size limit); 
finally, the intermediate state is a mixture of the two previous ones and it is found for 
e > e c and N < oo: both the consensus state and the coexistence one are locally stable, 
thus, the system is tri-stable, and the steady state depends on initial conditions. 

These results provide a nice illustration of how noise can effectively change the shape 
of the deterministic potential, allowing for noise-induced phenomena. Still, the presence 
of absorbing states - with the associated singularities in the steady state distribution 
- hinders true phase transitions to occur: the only possible steady state for any finite 
system is an absorbing one. Instead, in the infinite size limit, noise vanishes and the 
coexistence state becomes truly stable. 



The effect of quenched disorder in neutral theories 10 

3.3. Finite-N approximate solution: Introducing a "mutation" rate 

In order to regularize the singularities reported above and explore the possibility of 
phase transitions, we introduce a small "mutation" term. Mutation is defined as the 
process by which any randomly selected spin spontaneously inverts its state (regardless 
of its neighbors or external field) at some rate v > 0. The transition rates W become 

W* = (1 + e) [(1 - u){\ - x){\ + x - yj\ + |(| - x) 
W' x d = (1 - e) [(1 - v)x{\ - x + y)]+%x , , 

W* = (l + e)l{l-v){\-y){\+y-x)]+^\-y) l ] 

W« = (l- e )[(l-%(i-y + x)]+fy. 

Assuming v <C e <C 1 and keeping the leading orders in v and e, the Fokker-Planck 
equation for the global magnetization in presence of speciation becomes 

MM = -fy [AMMM\ + \% [BMMM] (is) 

with 

A v (<P) = -^<j)(l-<P 2 + 2v) 

(17) 
i3,(0) = [(l-e 2 )(l-0 2 ) + 2z,]/iV, 

respectively, and the associated stationary probability distribution function becomes 

1 (_K^ 

;i-e 2 )(l-0 2 ) + 2^ eXP V 2 1-e 2 
which, owing to u, does not have any singularity [51J. It is important to notice that if v 
is small enough, namely v ^C 2/(2 + N) (see Appendix), it does not affect significantly 
the dynamics, apart from removing the absorbing boundaries. 



V ^M - „ , 2Wl _^ 2 ^ 9 exp ---— ^ ] , (« 



We have used equation (18) to check our (approximate) results against numerical 
simulations of the complete (exact) dynamics (as obtained for a complete graph of 
iV spins by means of the Gillespie algorithm [S2]). Results for different values of e, 
are reported in Figure [3] which shows a rather good agreement with the theoretical 



predictions. Finally, from equation ( 18 ) it is easy to compute the value of e at which 



the second derivative of V v s computed in </> = changes sign in the small v limit: 

2TJV' < 19 > 

signaling a bifurcation from a unimodal distribution (coexistence) to a bimodal one. 
However, this is not a true thermodynamic phase transition; as expected, e c — > when 
iV — > oo. It is interesting to notice that the bifurcation from the absorbing to the active 
phase is completely 'noise-driven': any arbitrarily small amount of quenched-noise -i.e. 
any value of e > 0- leads to a stable active phase. Summing up: once absorbing states 
are perturbed with a non-vanishing mutation rate the only remaining stable-state in 
the thermodynamic limit is the active one, while for finite sizes there is a noise-induced 
transition. 
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Figure 3. Stationary probability distribution V s (<t>;e) in the case with mutation, 
v = 10~ 4 , for a finite complete network of size N = 200 and different values of 
the disorder strength, Dashed Lines: Curves for, e = 0.09 < e c ~ 0.1 (blue curve), 
e = 0.15 > e c (red curve), and e = 0.4 (green curve) as computed analytically from Eq. 



(18). Dots: Numerical results for simulations of the model (fully connected network 
and Gillespie algorithm) for the same parameters as above. 



4. Connection to the Generalized-voter class 



Al Hammal et al. [18] introduced a phenomenological Langevin equation aimed at 
capturing all the possible features of systems with two-symmetric absorbing states. 
Following the notation in [18] this equation reads 

<P(x,t) = V 2 4>(x,t) + {a<P - 60 3 )(1 - (j) 2 ) + y/l-^(x,t). (20) 

where a and b are constants, <f)(x,t) is a field whose dynamics is frozen if <p 2 = 1, and 
£(x, t) is a Gaussian white noise of zero mean and variance a. With the only requirement 
that b < this equation reproduces the critical behavior of the GV universality class. 
On the other hand, if b > then the equation resembles very much that of the Ising 
model (Model A), and indeed, in such a case one obtains that the GV transition is split 
into two different phase transitions: (i) the first one is Ising-like and corresponds to a 
breakdown of the up-down, Z 2 , symmetry; (ii) once the symmetry has been broken, by 
further changing parameters such that \fajb (fc [—1,1], the system eventually falls into 
the corresponding absorbing state (directed-percolation like phase transition). From 
this perspective, as pointed out in [18] the GV transition can be viewed as the merging 
of two different phenomena: the breaking down of a Z 2 symmetry and the falling down 
into an absorbing state, both of them occurring at the same transition point. 

If now we consider Eqs. ([7]) and (J8J) for the description of the disordered model 
considered here (without mutation), we can write it in the equivalent form of the 
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Langevin equation in the Ito prescription [H] 



6 2 



0(t) = --0(1 - 2 ) + y/(l-e*)(l-<P)r,(t) (21) 

where rj(t) is a Gaussian white noise with (rj(t)) = and (r)(t)r)(t')) = jrS(t — t'), 



which coincides with the O-dimensional version of eq. (20) (i.e. eq. (20) without spatial 

dependence) once the identifications b = 0, a = — e 2 /2 and a = \/^(l — e 2 ) are made. 
Therefore, at least at mean-field level, the VM with quenched random field closely 
resembles the GV dynamics (without quenched disorder). The main effect of random 
fields is to create a deterministic force which converts the state of coexistence (<f) = 0) 
into a stable one. Such a state is the only possible stationary state for infinitely large 
systems, while for finite-sizes there is a transition very similar to that of the GV class, 
without quenched disorder. Then, it is somewhat surprising that also this model is 
effectively well described by the same equations. 

It is also noteworthy that e = corresponds to the critical point in the 
thermodynamic limit and then, as e > 0, only the active phase exists: the absorbing 



phase of the AlHammars equation (20) is not accessible to the present model with 
quenched random fields. Thus, it is interesting to investigate what happens when 
quenched- disorder is introduced into a (non-linear) version of the VM, including an 
active and an absorbing phase, and a critical point separating the two of them. We 
tackle this problem in the next section. 

5. Nonlinear voter models, disorder, and spontaneous symmetry breaking 

So far we have seen that quenched disorder pushes a linear VM out of the criticality 
introducing an effective potential term that forces the system into an active symmetric 
phase. In this section we extend the analysis to the larger class of nonlinear voter models 
(NV) by mean field analysis and simulations of the model introduced in [53J in presence 
of a disordered environment. 

A general argument by Imry and Ma [13] predicts that, at equilibrium, quenched 
disorder prevents the spontaneous symmetry breaking of a discrete symmetry in D < 2 
and of a continuous symmetry in D < 4 [53] (a rigorous proof has been given later 
in [55]). However, in a remarkable recent paper Barghathi and Vojta [42J established 
that a one-dimensional model in the GV class violates the Imry-Ma result: owing to 
its non-equilibrium nature it can exhibit a spontaneous symmetry breaking even in low 
dimensions (D = 1 in their work). 

Thus, we want to study the effect of quenched random field disorder in the remaining 
case, i.e. when the non-linear voter model is such that it exhibits two separate transitions 
occurring separately. Does this scenario survive the presence of quenched disorder? 

We consider a particular nonlinear Voter Model introduced in [53] where, if x is 
the density of spins not aligned with Oi in its neighborhood, the flipping probability of 
a randomly selected spin is given by 

fi(x) = P(ai — > — o~i\x) oc xK(x), (22) 
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where K(x) is an arbitrary nonlinear function of x. We consider a cubic function of the 
form 

K{x) = ax 2 + ~b{-x~x 3 ) (23) 



that yields a stochastic equation of motion of the form of equation (20) with effective 
parameters a = a/4 and b = 6/16. If we choose b > 0, a > 0, this model has two 
(meta)stable fixed points at global magnetization 0* = ±^a/b [HJ [25] . Next, we 
include disorder by adding a quenched site-dependent term e, = — erjCj that locally 
breaks the up down symmetry favoring or disfavoring the single spin-flip. The quenched 
variables Tj are defined as in Section 2. 

5.1. Mean Field theory 

It can be conjectured that -in analogy with the linear case- for this nonlinear VM mean 
field in a random field, with parameters of nonlinearity a and b and disorder strength e, 



the presence of disorder leads to an equation of the form of Eq. (20), but with a new 
effective parameter a' = a — e 2 /2. Namely, we expect the system to be described by 

kt) = [(« - \)<t> - H 3 }^ - 2 ) + V(<f>, t). (24) 

From this equation we can compute the critical value of the external field intensity, 
e s c b , at which the disorder destroys the possibility of a spontaneous symmetry breaking, 
that is, ef = \/2a. Analogously, the position of the minimum of the potential (and 
therefore the maximum in the stationary probability distribution of the magnetization) 
is expected, when a, b > 0, at 



a-e 2 /2 , s 

0*^±Y f~. (25) 

These results turn out to be in good agreement with computer simulations in a fully 
connected network; in particular, in Figure [4] we show the dependence of the peaks 
position, 0*, on the intensity of the disorder and that it is independent of the size of the 
system. Observe that the up-down symmetry is spontaneously broken in this mean-field 
like case. 

5.2. Simulations in D = 2 

When D < 2, as discussed above, the phenomenology could be radically different if 
the Imry-Ma argument holds. In this case, the position of the peaks, p* = (0* + l)/2, 
depends on system size N. Indeed, in Figure [5] we show results for the position of the 
peaks as a function of 1/N in two different cases, with and without disorder. In the 
pure case (e = 0) the curve converges to some value between 0.25 and 0.3, implying 
that symmetry breaking is preserved in the thermodynamic limit. Instead, as soon as 
e ^ 0, the peak position tends to its symmetric-state value 1/2 as size is increased: 
the presence of disorder makes the symmetry breaking disappear for sufficiently large 
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Figure 4. Non-linear voter model running upon a fully connected network and 
quenched disorder of strength e: top: Stationary probability distributions P s {4>) plotted 
for e = 0.1 < €£* and various values of the system size, N. The position of the peak 
4>* does not depend on TV; the Ising up-down symmetry is broken, bottom: Position 
of the peak </>* for a — 0.0175/(1 + e) and b — 0.053/(1 + e) and varying e. The green 



dashed line represents the prediction Eq. ( 25 ) . In the inset are shown the probability 



distributions for some values of e: as e — > e c its two peaks become closer and closer, 
until eventually they merge at e c . 
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system sizes, in agreement with the argument of Imry and Ma for the equilibrium Ising 
symmetry breaking. 
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Figure 5. Two dimensional results for the non- linear voter model: position of the 
peaks as function of 1/N for the model with the same parameters as in figure El 
but in D = 2. The dashed line at p* = (0* + l)/2 = 0.5, to which the maxima 
of the distribution functions converge upon increasing system size, indicates that the 
symmetric active state is restored in the thermodynamic limit. 

Therefore, the non-linear voter model is a non-equilibrium one which-as predicted 
by the general theory of AlHammal et al.- can exhibit two different scenarios depending 
on parameter values: (i) either 2 separate phase transitions, Ising (for the up-down 
symmetry breaking) and DP (for the falling into the absorbing phase) or, instead, 
(ii) a unique GV transition (when up-down symmetry breaking and falling into the 
absorbing state occur simultaneously). In the first case, quenched disorder is a relevant 
perturbation, the Imry-Ma argument applies and no up-down symmetry breaking can 
occur in low dimensional systems (although it does occur at a mean- field level). Instead, 
in the second scenario, the GV is a genuine non-equilibrium phase transition, the Imry- 
Ma argument breaks down and the symmetry can be spontaneously broken even in 
low- dimensional cases. 

6. Conclusions 



Models with symmetric absorbing states mimicking neutral evolution -in its generalized 
sense- have been proposed in different areas of science, from physics and chemistry where 
they can model the behavior of kinetic reactions, to biology, genetics, and ecology where 
models like the voter model -also called Moran process- [56] or its dual representations 
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-as the coalescent- [57] have been successfully applied to understand (at least at a 
null-model level) different empirical observations. In all these models and applications, 
nevertheless, external and environmental "forces" are typically ignored. A qualitative 
change in the trustworthiness of the mathematical approach to these problems needs 
a stronger consideration of the fact that the particles/individuals/agents in the system 
do not act in the vacuum, but live in an environment that might both modify the 
interactions amongst individuals and interact themselves with the constituents of the 
system. As a first step towards a more detailed understanding of the effect of quenched 
disorder in such situations, in this paper we have studied how quenched random fields 
affect the dynamics in the voter model and in other similar non-linear models (in which 
the transition rates depend on the local state in a non-linear way). 

Our main findings are as follows: 

i) First, we investigated analytically the effects of a quenched random 'preference 
field' in the standard voter model dynamics, finding -in agreement with previous results- 
that the main effect of this type of noise is to generate a deterministic potential that 
stabilizes coexistence of opinions. We went further and found an effective stochastic 
description of the dynamics on a complete graph which is a function of the total number 
of spins N and the disorder strength e, and considered versions of the model allowing 
for spontaneous flip (mutations) from one state to its opposite. We report on disorder- 
induced transitions in finite systems, which disappear when the system size becomes 
infinite, i.e. in the thermodynamic limit; this occurs owing to the non-trivial interplay 
between noise and deterministic drift. 

ii) We have found a relation in terms of stochastic (Langevin) equations between a 
linear VM with disorder and the larger class of nonlinear voter models without disorder 
(the generalized voter Langevin equation, as proposed in [T8]). 

Hi) Finally, we have scrutinized the validity of the Imry-Ma argument in a family 
of non-equilibrium non-linear voter models. This was motivated by the recent finding 
by Barghathi and Vojta [12] that the Imry-MA argument can break down in non- 
equilibrium systems. In particular, if the corresponding pure model has a unique 
GV type of phase transition, then the corresponding disordered model which quenched 
random fields, can still exhibit a spontaneous breaking of the up-down symmetry even in 
low dimensions. This is a genuinely non-equilibrium feature, prohibited in equilibrium 
systems. 

On the other hand, as predicted by a general theory [18], some non-equilibrium 
non-linear voter models, do not undergo a unique phase transition but a sequence of 
two of them (one Ising like and one DP like). In this case, as we have shown here by 
means of computational studies and analytical arguments, the Imry-Ma criterion does 
hold and the up-down symmetry is not broken in low dimensional systems. 

The underlying reason for this, is that in the case studied by Barghathi and Vojta 
when the discrete up-down symmetry breaks down the system simultaneously falls into 
one of the two absorbing states; thus fluctuations are immediately extinguished. In the 
lack of stochasticity the Imry-Ma argument clearly breaks down as it rests in evaluating 
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the effects of eventual ffuctuations into ordered states. Instead in the other case studied 
here, the broken symmetric state to-be would be a fluctuating one and the corresponding 
phase transition an Ising one. Therefore, at sufficiently large scales, the system behaves 
as an equilibrium one, the Imry-Ma argument holds and no actual symmetry breaking 
occurs in low dimensions. 

We finally note that our description was limited for the sake of simplicity to only 
two species, but we do not expect qualitative changes for generalizations with more 
species [53J. 
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Appendix. Stationary distribution for the VM with mutation 

Consider a 2-species voter model upon the complete graph consisting of N nodes, 
without disorder, but with a mutation rate v > 0. Calling p the density of one of 
the two species, then the global magnetization is <fi = 2p — 1 and in analogy with 



equation (15) the birth-death rates describing the system's dynamics are 



W + = (i- v )(i-p) p + Z(i- p ) 

W. = (l- v )p(l-p) + -p, 

from which it is easy to obtain a Fokker-Planck equation for the probability distribution 
of the global magnetization, namely 

p(<M)= -94-|0P(0,t)}+ 

1 (27) 

+ ^dl{[(l-u)(l-^) + u]P(<p,t)}. 

Since the process is one-dimensional, it is straightforward [B] to write the equation for 
the stationary distribution of <fi, 

that can be imagined as derived from a potential V((f>;N,v). Given a system size N, 
the distribution, initially peaked at the boundaries = ±1, starts to have a peak at 
= when 

' (29) 



2 + iV' 
that can be easily verified computing the second derivative <9?Poo(0; N, u)\^ = o. 
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Thus, for finite-size systems, when v <C 2/(2 + N) the introduction of a mutation 
rate has the effect of removing the absorbing points of the process and it justifies the 



assumption in section 3J3 It is worth noting that the Langevin equation of this process, 
with mutation but without disorder, is very close to equation ^ if v = e 2 . This 
backs the intuition that the quenched external field acts effectively as a diffuse source of 
particles (individuals) in the same way as a mutation term does, but with the substantial 
difference that while mutation eliminates the absorbing states, this does not happen in 
the case of the disordered system. 
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